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ABSTRACT 


Several  computer  programs  for  the  analysis  of  longitudinal  wave  cuts 
which  were  developed  at  the  University  of  Michigan  under  the  direction  of 
S.D.  S harms  have  been  converted  at  the  David  W.  Taylor  Naval  Ship  Research 
and  Development  Center  (DTNSRDC)  for  use  on  the  CDC  6000  Computer  Series. 
The  computer  programs  that  are  documented  In  this  report  perform  the  fol- 
lowing tasks:  convert  digital  wave  cut  data  for  use  In  a wave  analysis 
program  which  determines  free-wave  spectra  for  a number  of  transverse  wave 
numbers  and  calculates  the  wave  resistance;  plot  the  free  wave  spectra; 
determine  and  plot  the  contours  of  bulb  Influence  factors.  The  programs 
have  been  tested  at  DTNSRDC  on  a sample  of  wave  cut  data  measured  at  the 
University  of  Michigan. 
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ADMINISTRATIVE  INFORMATION 


The  work  reported  herein  w»3  funded  under  NAVSSC  Project  Order  60IA2. 

The  work  was  performed  under  David  W.  Taylor  Naval  Ship  Research  and 
Development  Center  Work  Unit  1524-599. 

INTRODUCTION 

The  Naval  Ship  Engineering  Center  (NAVSEC)  requested  that  the  David  W. 
Taylor  Naval  Ship  Research  and  Development  Center  (DTNSRDC)  convert  several 
computer  programs  for  use  on  the  CDC  6000  computer  serias  and  develop  the 
necessary  documentation  to  use  and  maintain  these  programs.  These  computer 
programs  were  developed  at  The  University  of  Michigan  under  the  direction  of 
S.  D.  Shanna. 

Collectively,  these  programs  are  used  for  the  analysis  of  longitudinal 
wave  cuts.  Individually,  the  programs  perform  the  following  tasks:  the 
first  program  NONDIM  converts  digital  wave  cut  information  for  use  by  the  wave 
analysis  program;  the  Wave  Analysis  Program  (WAVECT)  performs  wave  analysis  to 
determine  the  free-wave  spectra  for  a number  of  transverse  wave  numbers*  with 
an  optional  truncation  correction  for  height  data  and  also  calculates  the 
wave  resistance;  the  Free-Wave  Spectra  Plotting  Program  (SPCTRA)  is  used  to  plot 
the  free  wave  spectra  derived  by  the  WAVECT  program;  the  Bulb  Con*-cur  Printer 
Program  (BUBOPT)  uses  the  free-wave  spectra  provided  by  the  WAVECT  program  to 
determine  the  contours  of  the  bulb  influence  factor  for  use  by  the  Bulb  Contour 
Plotter  program  (BUBPLT) ; the  BUBPLT  program  plots  the  contours  of  the  constant 
bulb  influence  factor  (ETA). 

The  above  programs  have  been  converted  and  tested  on  the  CDC  6000  Computer 
Series  at  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center.  The 
testing  used  data  reported  in  Sharma  and  Naegle  (1970).* 

*The  nondimensional  transverse  wave  number  is  defined  as:  u = Sec  9 Tan  9, 
where  the  angle  9 defines  the  direction  of  wave  propagation. 

^References  are  listed  on  page  76. 
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BACKGROUND 

For  many  years  there  has  been  continuing  interest  in  the  reduction  of 

resistance  into  components  attributable  to  specific  causes  such  as  friction, 

wind,  and  wave  breaking.  In  recent  years  the  development  of  theory  and 

experimental  techniques  have  made  it  possible  to  assess  more  realistically 

the  contribution  of  wave  making  phenomena  to  the  total  resistance  of  ships. 

Specifically,  the  use  of  the  longitudinal  wave  cut  technique  with  theory  to 

determine  the  wave  making  characteristics  of  ship  hulls  from  model  experiments 

2 

has  been  advanced  by  a number  of  investigations. 

In  1968  S.  D.  Sharma  directed  a project  at  the  University  of  Michigan 
with  the  objective  of  optimizing  bulbous  bows  for  ships.  The  approach  was 
to  use  experimental  methods  to  determine  the  profiles  of  wave  patterns  for 
specific  hull  and  bulb  configurations  and  then  use  theoretical  methods  to 
predict  the  effects  of  changes  in  bulb  size  and  location  on  the  wave  patterns 
and,  consequently,  on  the  component  of  resistance  due  to  wave  making.  Wave 
profile  measurements  were  obtained  for  three  bulbs  fitted  to  a common  hull 
model.  The  investigators  concluded  that  the  technique  was  promising  as  a 
design  tool  for  optimizing  bulb  size  and  location  using  a minimum  of  model 
experiments. 

Since  the  methods  developed  by  Sharma  appear  to  have  promise  in 

optimizing  hull  designs,  the  programs  have  been  adapted  for  use  by  the 

U.S.  Navy  as  described  herein.  These  programs  have  also  been  used  in 

3 

modified  form  by  other  investigators. 

COMPUTER  PROGRAM  DESCRIPTIONS  (USER'S  GUIDE) 

This  section  of  the  report  provides  a user's  guide  for  the  five  computer 
programs  which  collectively  will  be  used  for  the  analysis  of  longitudinal  wave 
cuts.  The  program  operating  instructions  are  for  the  CDC  6000  computer  series. 
The  descriptions  provided  summarize  the  purpose,  input  requirements /format  and 
output  for  each  of  the  computer  programs. 

NOND IMEN S IONAL  WAVE  CUT  PROGRAM  (NONDIM) 

The  program  NONDIM  converts  digital  wave  cut  information  from  analog 
units  (i.e.,  millivolts/inch)  to  feet,  and  then  nondimensionalizes  the  wave 
cut  information  for  use  in  the  Wave  Analysis  Program  (WAVECT) . 
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The  input  variables  for  program  NONDIM  are  defined  as  follows: 

LABEL  - A vector  containing  the  80  character  title 

MMU  - A control  character  indicating  whether  height  data  or 
slope  data  is  to  be  read  in: 

MMU  j<  0;  Height  data 

MMU  > 0;  Slope  data 

MAX  - The  number  of  points  in  the  wave  profile,  ZETA  - MAX 
must  be  less  than  or  equal  to  1000 

MEAN  - The  number  of  points  in  the  wave  profile  which  should  be 

averaged  and  subtracted  from  all  of  the  points  of  the  wave 
profile  in  order  to  obtain  a zero  line 

DT  - The  time  interval  between  successive  input  points  in  the 
profile  in  seconds 

YFT  - The  transverse  distance  from  the  center  line  of  the  model 
to  the  center  wire  of  the  wave  probe  in  feet 

VFPS  - The  model  speed  in  feet  per  second 

ZSCA  - Scale  factor  for  use  in  converting  the  input  data  into 

either  inches  or  slope  - ZSCA  is  the  number  of  input  units 
per  inch  of  wave  height,  or  per  radian  of  slope  as  indicated 
by  MMU. 

XSHUT  - Distance  from  the  origin  on  the  model  to  the  trip  in  feet, 
with  positive  measured  forward  - See  Figure  1 

XLAG  - Distance  from  center  wire  of  the  wave  probe  to  the  center 
of  the  photo  cell  in  feet  with  positive  upstream  - See 
Figure  1 

JMARK  - The  number  of  time  intervals  (fractions  allowed)  from 

the  trip  to  where  the  input  begins,  with  positive  measured 
forward  - See  Figure  2 

ZETA  - A vector  containing  the  MAX  input  points  for  the 
longitudinal  wave  profile. 
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The  input  cards  for  the  NONDIM  program  are  read  from  logical  device 
Number  5 and  should  be  prepared  in  the  following  manner: 


CARD 

VARIABLES 

FORMAT 

1 

LABEL 

8A10 

2 

MMU,  MAX,  MEAN 

315 

3 

DT,  YFT,  VFPS 

3F10.3 

4 

ZSCA,  XSHUT,  XLAG, 

JMARK 

4F10.3 

5 

IFMT 

8A10 

6 

ZETA  (As  many  cards 
as  necessary) 

IFMT 

Below  is  a list  of  input  values  used  to  complete  one  successful  run 
of  the  NONDIM  program.  Because  of  the  size  of  ZETA,  the  specific  input 
values  for  ZETA  are  omitted  from  the  list  but  are  provided  in  Table  1. 


LABEL  - Model  1094  (CV)  RUN  NO.  4 AUG.  6,  1978  V = 5.360  FT/SEC 

MMU  = 0 

MAX  * 451 

MEAN  - 50 

DT  » 0.03 

YFT  = 4.125 

VFPS  - 5.360 

ZSCA  - -3.972 

XSHUT  - 3.0833 

XLAG  - 1.1667 

JMARK  - 118.25 

IFMT  - (15F5.2,  5X) 

Table  1 provides  a listing  of  the  above  sample  input  in  the  format 
required  by  the  progran  NONDIM. 
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The  output  from  the  NONDIM  program  consists  of  two  types;  formatted 
printer  output  and  unformatted  file  output.  The  formatted  output  consists 
of  the  longitudinal  step  size  (DX)  which  is  the  distance  between  input 
points  in  the  wave  profiles  (DX  should  be  negative) ; the  distance  of  the 
first  input  point  forward  from  the  origin  on  the  model  (XONE) ; the  distance 
from  the  centerline  of  the  model  to  the  center  of  the  wave  wire  (Y) ; and 
the  number  of  points  in  the  wave  profile  (MAX).  This  is  followed  by  the 
points  of  the  wave  profile  (ZETA) , which  have  been  adjusted  to  zero  by 
averaging  a number  of  points  from  the  wave  profile  (this  number  is  determined 
by  the  input  variable  MEAN) . Both  dimensional  and  nondimens ional  forms 
are  printed.  Table  2 provides  a sample  listing  of  the  dimensional  output; 
Table  3 provides  a sample  listing  of  the  nondimens ional  output.  The  non- 
dimensionalized  data  is  also  written  unformatted  on  a file  for  use  by  the 
Wave  Analysis  Program  (WAVECT)  and  the  Free-Wave  Spectra  Plotting  Program 
(SPCTRA).  The  length  information  is  nondlmensionalized  by  the  fundamental 
wave-number  [g/u2]. 

The  NONDIM  program  uses  device  numbers  4,  5,  and  6 for  its  input  and 
output.  Device  number  5 should  be  assigned  to  the  card  reader  for  the 
input,  device  number  6 should  be  assigned  to  the  line  printer  for  the 
formatted  output,  and  device  number  4 should  be  assigend  to  a permanent 
file  (disk  or  magnetic  tape)  which  needs  to  be  requested  and  catalogued. 


Sample  control  cards  used  to  run  program  NONDIM  on  the  CDC  6000  Series 
Computers  are  as  follows: 

CHRENON , CM60000 , T100 , P3 . 

CHARGE , CHRE , XXXXXXXXXX . 

REQUEST, TAPE4,*PF. 

FTN(T) 


CATALOG , TAPE4 , PERMFILE4 , ID-CHRE . 
7/8/9  END  OF  RECORD  CARD 
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7/8/9  END  OF  RECORD  CARD 
Data  Cards 

6/ 7/8/9  END  OF  FILE  CARD 

Note:  PERMFILE4  is  the  name  of  the  permanent  file  (disk)  on  which  the 
unformatted  output  is  to  be  stored. 

WAVE  ANALYSIS  PROGRAM  (WAVECT) 

The  program  WAVECT  determines  the  free-wave  spectra,  wave  resistance, 

and  side  force  on  a body  from  a wave  profile  measured  along  a line  parallel 

to  the  path  of  the  model  (longitudinal  wave  cut) . The  method  used  is  that 

2 

developed  by  Sharma.  WAVECT  uses  the  nondimensionalized  data  from  the 
program  NONDIM  and  does  the  actual  wave  analysis  with  an  optional  truncation 
correction  which  can  be  applied  to  height  data. 

The  input  variables  for  the  WAVECT  program  are  defined  as  follows: 


MU  - Control  character  indicating  whether  height  or  slope  data  is 
to  be  read  and  whether  or  not  the  truncation  correction  is  to 
be  made  in  the  case  of  height  data: 

MU  < 0;  height  data  without  truncation  correction 

MU  = 0;  height  data  with  truncation  correction 

MU  > 0;  slope  data 

N - The  number  of  elementary  waves  to  be  analyzed  during  the  wave 
aralysis. 

DU  - The  size  of  the  increment  for  the  transverse  wave  number  which 
is  used  to  determine  the  N elementary  waves  which  are  to  be 
analyzed 

M - The  number  of  points  which  are  to  be  read  from  the  wave  profile 

MP  - The  number  of  the  points  on  the  wave  profile  at  which  the  analysis 
for  determining  the  asymptotic  behavior  of  the  wave  profile  should 
begin  - This  is  used  in  making  the  truncation  correction  and  is 
not  read  unless  MU  ■ 0. 


1 


/ 


C3  -A  nondimenslonal  distance  indicating  where,  in  relation  to  the 
origin  on  the  model,  is  the  origin  for  use  in  determining  the 
beginning  of  the  asymptotic  behavior  of  the  wave  - This 
variable  is  not  read  unless  MU  * 0. 

LABEL  - A vector  which  contains  the  80  character  title 

DX  - Distance  between  the  input  points  in  the  wave  profile 

XONE  - The  distance  of  the  first  input  point  forward  of  the  origin 
on  the  model 

Y - The  distance  from  the  centerline  of  the  model  to  the  center 

of  the  wave  wire 

MAX  - The  number  of  points  in  the  wave  profile 

MMU  - A control  character  indicating  whether  height  or  slope  data 
is  to  be  read  in: 

MMU  0;  Height  data 

MMU  > 0;  Slope  data 

ZETA  - A vector  containing  the  MAX  input  points  for  the  longitudinal 
wave  profile. 


The  input  variables  LABEL,  DX,  XONE,  Y,  MAX,  MMU,  and  ZETA  are  the 
unformatted  output  of  Program  NONDIM  and  are  read  by  Program  WAVECT  from 
device  number  4.  The  input  variables  MU,  N,  DU,  M,  MP,  and  C3  are  user 
Inputs  read  by  WAVECT  from  device  number  5 and  should  be  prepared  as 
follows: 


CARD 

1 

2 

3 (if  MU  j*  0) 
3 (if  MU  - 0) 


VARIABLES 

MU 

N,  DU 
M 

M,  MP,  C3 


FORMAT 

15 

15,  F10.3 
15 

215,  F10.3 


Listed  below  are  sample  input  values  used  to  complete  a successful  run 
of  the  WAVECT  program. 
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MU  - 0 


i 


! 


N - 100 

DU  - 0.1 

M - 451 

MP  - 416 

C3  - 0.0 

Table  4 provides  a listing  of  the  above  sample  Input  In  the  format 
required  by  program  WAVECT. 

The  output  of  the  Wave  Analysis  Program  consists  of  a number  of  error 
statements,  one  of  two  forms  of  printed  output  on  device  number  6,  and 
the  free-wave  spectra  written  on  device  number  3 which  is  a disk  file, 
for  use  by  the  Bulb  Contour  Printer  program.  The  first  form  of  the  printed 
output  is  that  generated  for  either  slope  data  or  height  data  without  the 
truncation  correction,  the  second  output  form  is  generated  when  the  cal- 
culations on  height  data  are  made  with  the  truncation  correction. 

The  error  statements  are  printed  to  indicate  that  there  is  an  error 
in  the  input,  or  that  the  user  input  is  inconsistent  with  the  input  which 
is  read  from  the  Nondimensionalization  Program  (NONDIM)  output.  Error 
statements  are  output  by  program  WAVECT  if  any  of  the  following  conditions 
exist: 

1)  If  either  N or  DU  is  less  than  or  equal  to  zero 

2)  If  MU  indicates  a different  type  of  data  than  MMU  from  the 
Nondimensionalization  program  (i.g.,  if  MU  = +1  and  MMU  = 0, 
the  Wave  Analysis  Program  expects  slope  data  while  the  Non- 
dimensionalization program  is  giving  height  data) 

3)  If  M is  not  between  one  and  MAX 

4)  If  MP  is  not  between  one  and  M (this  applies  only  for  height  data 
with  truncation  correction) . 

Execution  of  the  program  is  terminated  after  any  of  the  above  errors  have 
been  detected. 
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If  no  errors  have  been  found,  the  program  then  proceeds  to  calculate  the 
free-wave  spectra,  transverse  force,  and  wave  resistance,  and  to  print  the 
output.  The  printer  output  on  device  number  6 consists  of  the  title,  the 
number  of  elementary  waves  to  be  analyzed,  the  transverse  wave  number  step 
size,  the  longitudinal  distance  between  the  points  in  the  wave  profile, 
and  the  number  of  points  In  the  wave  profile.  If  the  truncation  correction 
Is  to  be  made  (height  data  only)  the  point  In  the  wave  profile  where  the 
asymptotic  behavior  of  the  wave  is  assumed  to  begin  Is  printed  along  with 
the  three  constants  describing  the  asymptotic  behavior  of  the  wave.  This 
Is  then  followed  by  the  definition  of  all  terms  given  In  the  heading  of  the 
output  columns.  These  headings  and  their  definitions  for  the  output  with 
slope  data  and  height  data  without  the  truncation  correction  are  as  follows: 

NU  - The  number  of  the  elementary  wave  being  analyzed 

UNU  - The  transverse  wave  number  - Sec  9 • Tan  0 ^ 

SNU  - The  longitudinal  wave  number  - Sec  0 


/ 


TNU  - Tan  9 

SNU*  - The  sine  component  of  the  Fourier  transform  of  the  wave  profile 
CNU*  - The  cosine  component  of  the  Fourier  transform  of  the  wave  profile 

F - Hie  sine  component  of  the  free-wave  spectra  calculated 

G - The  cosine  component  of  the  free-wave  spectra  calculated 
E - The  square  root  of  the  sum  of  the  squares  of  F and  G 

T - The  nondlmenslonal  transverse  force •♦ 

R - The  nondlmenslonal  wave  resistance. 


♦♦The  wave  pattern  on  each  side  of  the  ship  exerts  a transverse  force  on  the 
ship.  If  the  ship  Is  symmetric,  It  generates  a symmetric  wave  pattern,  and 
therefore  antisynmetrlc  transverse  forces.  If,  however,  the  ship  Is  asym- 
metric, It  generates  an  asynmetrlc  wave  pattern,  which  exerts  an  asymmetric 
transverse  force  on  the  ship.  The  net  transverse  force  on  an  asymmetric  ship 
due  to  wave  making  can  be  determined  by  taking  the  difference  between  the 
side  forces  calculated  from  the  wave  pattern  measured  on  each  side  of  the 
ship.  The  longitudinal  force  on  an  asymmetric  ship  Is  the  average  of  the 
wave  resistance  computed  from  the  wave  pattern  on  each  side  of  the  ship. 


v , 
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The  variables  T and  R,  printed  with  each  step,  are  the  integrals  of  the 
free-wave  spectra  through  that  elementary  number. 


When  the  wave  analysis  is  applied  to  height  data  with  truncation 
correction,  the  output  on  device  number  6 is  as  follows: 

NU  - The  number  of  the  elementary  wave  being  analyzed 

UNU  - The  transverse  wave  number  - Sec  9 • Tan  0 

SNU  - The  longitudinal  wave  number  - Sec  9 

TNU  - Tan  0 

SNU*  - The  sine  component  of  the  Fourier  transform  applied  to  the 
wave  profile  without  the  truncation  correction 

SNU*T  - The  sine  component  of  the  Fourier  transform  of  the  wave  profile 
with  the  truncation  correction 

CNU*  - The  cosine  component  of  the  Fourier  transform  of  the  wave 
profile  without  the  truncation  correction 

CNU*T  - The  cosine  component  of  the  Fourier  transform  of  the  wave 
profile  with  the  truncation  correction 

F - The  sine  component  of  the  free-wave  spectra  calculated  with 
the  truncation  correction 

G - Hie  cosine  component  of  the  free-wave  spectra  calculated  with 
the  truncation  correction 

E - The  square  root  of  the  sum  of  the  squares  of  F and  G 

T - The  nondimens ional  transverse  force  calculated  without  the 

truncation  correction* 

TT  - The  nondlmensional  transverse  force  calculated  with  the 
truncation  correction* 

R - The  nondlmensional  wave  resistance  calculated  without  the 
truncation  correction 

RT  - The  nondlmensional  wave  resistance  calculated  using  the 
truncation  correction. 

The  wave  resistance  and  side  force  are  nondimensionalized  by  multiplication 
by  l/Pg"2V6. 

A sample  output  listing  showing  the  output  of  program  WAVECT  is 
provided  in  Table  5. 


*See  the  footnote  on  page  16. 
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thf.  nonoimensipnal  transverse  force  without  truncation  correction 


Table  5 - Sample  Output  for  Program  WAVECT  (Cont'd) 


002S  .0026 


Table  5 - Sample  Output  for  Program  WAVECT  (Cont’d) 


The  input  and  output  of  the  WAVECT  program  is  done  using  devices 
number  3,  4,  5,  and  6.  Device  5 should  be  assigned  to  the  card  reader 
for  the  user  input,  device  6 should  be  assigned  to  the  line  printer  for 
the  printed  output,  device  4 should  be  assigned  to  permanent  file 
containing  the  unformatted  nondlmenslonal  output  of  program  NONDIM  to 
be  used  as  input,  and  device  3 should  be  assigned  to  a permanent  file 
which  is  requested  and  catalogued  for  the  free-wave  spectra  output  of 
WAVECT.  The  file  containing  the  free-wave  spectra  output  is  used  by  the 
Bulb  Contour  Printer  Program  (BUBOPT)  and  the  Bulb  Contour  Plotter  Program 
(BUBPLT) . 

Sample  control  cards  used  to  run  program  WAVECT  on  the  CDC  6000 
Series  Computers  are  shown  below: 

CHREFWS , CM46000 , T100 , P3 . 

CHARGE , CHRE , XXXXXXXXXX . 

REQUEST, TAPE3,*PF 

ATTACH, TAPE4 ,PERMFILE4 , ID=CHRE . 

FTN(T) 

LGO. 

CATALOG , TAPE3 , PERMFILE3 , ID=CHRE . 

7/8/9  END  OF  RECORD  CARD 

Source  Deck 

7/8/9  END  OF  RECORD  CARD 

Data  Cards 

6/ 7/8/9  END  OF  FILE  CARD 

Note:  PERMFILE4  is  the  name  of  the  permanent  file  which  contains  the 

unformatted  nondlmenslonal  output  of  program  NONDIM,  PERMFILE3  is 
the  name  of  the  permanent  file  on  which  the  free-wave  spectra 


FREE  WAVE  SPECTRA  PLOTTING  PROGRAM  (SPCTRA) 


V 


k 


The  SPCTRA  program  plots  the  sine  component  of  the  free-wave  spectra 
(F) , the  cosine  component  of  the  free-wave  spectra  (G) , and  the  amplitude 
(E). 

The  input  variables  for  the  SPCTRA  program  are  defined  as  follows: 


N - The  number  of  free-wave  spectra. 

DU  - The  transverse  wave  number  step  size. 


LABEL  - An  80  character  title  for  the  output. 


F - The  sine  component  of  the  free-wave  spectra. 

G - The  cosine  component  of  the  free-wave  spectra. 


The  user  inputs  (N,  DU,  and  LABEL)  are  read  from  device  number  5 and 


should  be  prepared  as  follows: 

CARD 

VARIABLES 

FORMAT 

1 

N,  DU 

15,  F10.5 

2 

LABEL 

8A10 

The  variables  F and  G are  outputs  of  the  Wave  Analysis  Program  (WAVECT) 
and  are  read  by  SPCTRA  from  device  number  4. 

The  values  of  user  input  shown  below  were  used  to  complete  a 
successful  run  on  the  CDC  6000  Series  Computer. 

N - 100 

DU  =0.1 

LABEL  = MODEL  1094  (CV)  RUN  NO.  4 AUG.  6,  1968  V = 5.360  FT/SEC 

Table  6 provides  a listing  of  the  above  sample  input  in  the  format 
required  by  the  program  SPCTRA. 


I 


I 
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The  output  for  the  SPCTRA  program  consists  of  plots  of  the  sine 
component  of  the  free-wave  spectra  (F)  versus  speed  (S)  and  versus 
heading  (U) , the  cosine  component  of  the  free-wave  spectra  (G)  versus 
S and  versus  U,  and  for  the  amplitude  (E  = ^F2  + G2)  versus  S and  versus  U. 
F,  G,  and  E are  plotted  along  the  vertical  axis  while  U and  S vary  along 
the  horizontal  axis.  The  plot  file  is  automatically  written  on  a magnetic 
tape  which  has  to  be  requested  in  the  control  cards.  The  tape  is  then 
mounted  on  the  CalComp  Plotter  which  does  the  actual  plotting.  Figure  3 
is  a sample  of  the  plotted  output  of  program  SPCTRA  using  the  CalComp 
plotter. 

The  input  and  output  of  the  SPCTRA  program  is  done  via  devices  numbers 
4,  5,  6,  and  7.  Device  number  5 should  be  assigned  to  a card  reader  for 
user  inputs.  The  printed  output  is  printed  on  logical  device  number  6. 
Device  number  4 should  be  assigned  to  the  permanent  file  written  by  the 
Wave  Analysis  Program  (WAVECT)  to  be  used  as  input  and  device  number  7 
should  be  assigned  to  a 7-track  magnetic  tape  on  which  the  commands  are 
written  which  drive  the  CalComp  plotter. 

Sample  control  cards  used  to  run  SPCTRA  on  the  CDC  6000  Series 
Computers  are: 

CHREFWS , CM50000 ,MT1 , T100 , P3 . 

CHARGE , CHRE , XXXXXXXXXX . 

FTN(T) 

ATTACH , TAPE4PERMFILE4 , ID=CHRE . 

VSN,TAPE7=RAE01=SL0T32. 

REQUEST ,TAPE7 ,HI .RING.  (SLOT 32 .RAE01) 

ATTACH, CALC9 36. 

LDSET, (LIB*CALC936) 

LGO. 

RETURN, TAPE7. 

7/8/9  END  OF  RECORD  CARD 

Source  Deck 
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AUG.  6,  1968  V-5.360  FT/ SEC 


ile  Output  for  Program  SPCTRA:  Sine  Component  (F) , Cosine  (G) 
Amplitude  (El  versus  Transverse  Wave  Number  (Sec  0 TAN  0) 


MODEL  1094  (CV)  RUN  NO.  4 AUG.  6,  1968  V=5.360  FT/SEC 


7/8/9  END  OF  RECORD  CARD 


l ! 


I 

l 


Data  Cards 

6/7/ 8/ 9 END  OF  FILE  CARD 

Note:  PERMFILE4  is  the  name  of  the  permanent  file  on  which  the  free-wave 
spectra  output  of  program  WAVECT  resides,  RAE01  is  the  name  of  the 
7-track  magnetic  tape  being  used  which  must  be  provided  by  the  user. 

BULB  CONTOUR  PRINTER  PROGRAM  (BUBOPT) 

The  program  BUBOPT  determines  the  coarse  contours  of  the  bulb  influence 
factor.  This  program  is  primarily  Intended  for  determining  the  range  of 
variables  to  use  as  input  to  the  BUBPLT  (Plotter)  program. 


The  input  variables  used  in  the  BUBOPT  program  are  defined  as  follows: 
LABEL  - An  80  character  title  for  the  output 
DU  - The  transverse  wave  number  step  size 
L - The  nondimens ional  length 

LAMBDA  - A control  character  indicating  the  type  of  data  to  be  read: 

LAMBDA  = 0;  hull-bulb  spectra  is  given  and  the  bulb  spectra 
must  be  separated  from  this 

LAMBDA  ^ 0;  bulb  spectra  given 

N - The  number  of  free-wave  spectra 

NCONT  - A control  character  controlling  the  calculation  of  contour 
points : 

NCONT  = 0;  calculate  and  print  the  points  for  plotting 
contours  of  NU  • constant 

NCONT  i4  0;  do  not  calculate  points  for  the  contours 
QMIN  - Starting  location  of  bulb  as  a fraction  of  ship  length* 

QMAX  - Ending  location  of  bulb  as  a fraction  of  ship  length* 

DQ  - The  step  size  for  Q 

*Positive  locations  mean  the  new  bulb  is  forward  of  the  location  of  the 
current  bulb,  and  negative  locations  are  aft. 


t 


ETAMIN  - The  minimum  eta  for  which  contours  are  wanted  (must  be 
zero  or  larger)  “ Eta  is  defined  as 


w 


OPT 


w. 


BH 


ETAMAX  - The  maximum  eta  for  which  contours  are  wanted 
DETA  - The  distance  between  contours 


FH 


- Sine  component  of  the  free-wave  spectra  without  bulb  output 
by  program  WAVECT 


GH 


- Cosine  component  of  the  free-wave  spectra  without  bulb 
output  by  program  WAVECT 


FBO 


- Sine  component  of  the  free-wave  spectra  with  bulb  output  by 
program  WAVECT 


GBO 


- Cosine  component  of  the  free-wave  spectra  with  bulb 
output  by  program  WAVECT. 


The  user  input  variables  LABEL,  DU DETA  for  the  BUBOPT  program 

are  read  from  device  5 and  should  be  prepared  as  follows: 


CARD 


VARIABLES 


FORMAT 


LABEL 


DU,  L,  LAMBDA, 
N,  NCONT 


8A10 

2F10.5,  315 


QMIN,  QMAX,  DQ, 
ETAMIN,  ETAMAX, 
DETA 


6F10.5 


Values  of  user  input  data  that  were  used  to  perform  a sample  run  of 
BUBOPT  on  the  CDC  6000  Series  Computer  are  provided  below: 


LABEL  - FROM  MODELS  1094CV  AND  1094-B2  AT  V = 5.360  FT/SEC  (AUG.  1968) 
DU  - 0.1 

L - 14.01 

LAMBDA  - 0 

N - 100 
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NCONT 

- 0 

QMIN 

- -.08 

QMAX 

- 0.08 

DQ 

- 0.01 

ETAMIN 

- 0.6 

ETAMAX 

- 2.0 

DETA 

- 0.1 

Table  7 provides  a listing  of  the  above  sample  input  in  the  format 
required  by  program  BUBOPT. 

The  input  variables  FH  and  GH  are  read  from  device  number  3 and  should 
be  the  free-wave  spectra  output  of  program  WAVECT  without  a bulb.  Input 
variables  FBO  and  GBO  are  read  from  device  number  4 and  should  be  the 
free-wave  spectra  output  of  program  WAVECT  with  a bulb. 


i 


/ 


i 


The  output  of  the  BUBOPT  (Printer)  program  consists  of  the  following: 

Q - The  longitudinal  location  of  the  bulb 
P(MIN)  - P location  of  ETA(MIN)  for  a given  Q 
A,  B,  C - For  a given  Q;  ETA  is  quadratic  in  P 
ETA  - AP2  + BP  + C 
P(l)  - The  bulb  volume  for  eta  * 1 

For  eta  = constant,  a given  Q gives  PI  and  P2  for  plotting  contours. 

It  obtains  PI  and  P2  from  the  solution  to  AP2  + BP  + (C  - ETA) 

= 0,  for  eta  = constant.  PI  and  P2  = -B+Vb^  - 4A  (C  - ETA) 

2A 

P - Bulb  volume  - current  bulb  is  (P=l)  - As  fraction  of  bulb 
volume 

Table  8 is  a sample  of  the  output  and  output  format  that  is  generated 
by  program  BUBOPT. 

The  input  and  output  of  program  BUBOPT  is  done  via  devices  3,  4,  5, 
and  6.  Device  number  5 should  be  assigned  to  the  card  reader  for  the  user 
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29 


Table  7 - Sample  Input  for  Program  BUBOPT 


Program  BUBOPT 


Table  8 - Sample  Output  for  Program  BUBOPT  (Cont’d) 


i 


! 


: 
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’ — ...'m 

I 

inputs,  device  6 should  be  assigned  to  the  line  printer  for  the  printed 
output,  device  3 should  be  assigned  to  the  permanent  file  containing  the 
free-wave  spectra  output  of  Program  WAVECT  without  a bulb,  and  device  4 
should  be  assigned  to  the  permanent  file  containing  the  free-wave  spectra 
output  of  program  WAVECT  with  a bulb. 

Sample  control  cards  to  run  the  program  BUBOPT  on  the  CDC  6000  Series 
Computers  are  shown  below: 

CHREOPT , CM60000 , T120 , P 3 . 

CHARGE , CHRE , XXXXXXXXXX . 

ATTACH ,TAPE3 .PERMFILE3 , ID=CHRE . 

ATTACH , TAPE4 , PERMFILE4 , ID=CHRE . 

FTN. 

LGO. 

7/8/9  END  OF  RECORD  CARD 
Program  Cards 
7/8/9  END  OF  RECORD  CARD 
Data  Cards 

6/ 7/8/9  END  OF  FILE  CARD 

Note:  PERMFILE3  and  PERMFILE4  are,  respectively,  the  names  of  the  permanent 
files  containing  the  free-wave  spectra  output  of  program  WAVECT  with 
and  without  a bulb. 

BULB  CONTOUR  PLOTTER  PROGRAM  (BUBPLT) 

The  program  BUBPLT  determines  and  plots  the  contours  of  the  constant 
bulb  influence  factor  (ETA)  in  the  P-Q  space. 

The  input  variables  for  the  BUBPLT  program  are  defined  as  follows: 
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LABEL  - A title  for  the  output,  up  to  80  characters 
NCHAR  - The  number  of  characters  in  the  title 

N - The  number  of  free-wave  spectra 

LAMBDA  - A control  character  indicating  the  type  of  data  to  be  read  in 

LAMBDA  = 0;  the  hull-bulb  spectra  is  given,  and  the  bulb 
spectra  must  be  separated  from  this 

LAMBDA  t 0;  bulb  spectra  given 

DU  - The  transverse  wave  number  step  size 

L - The  nondimensional  length 

QMIN  - Starting  location  of  bulb  as  a fraction  of  ship  length 
QMAX  - Ending  location  of  bulb  as  a fraction  of  ship  length 
DQ1  - The  step  size  for  Q 

ETAMIN  - The  minimum  eta  for  which  contours  are  wanted 

ETAMAX  - The  maximum  eta  for  which  contours  are  wanted 

DETA1  - The  distance  between  contours 

DETA2  - The  distance  between  contours 

FH  - Sine  component  of  free-wave  spectra  without  bulb  output  by 
program  WAVECT 

GH  - Cosine  component  of  free-wave  spectra  without  bulb  output 
by  program  WAVECT 

FBO  - Sine  component  of  free-wave  spectra  with  bulb  output  by 
program  WAVECT 

GBO  - Cosine  component  of  free-wave  spectra  with  bulb  output 
by  program  WAVECT 

The  user  input  variables  LABEL,  NCHAR,  . . . , DETA2  for  the  BUBPLT 
program  are  read  from  device  5 and  should  be  prepared  as  shown  below: 
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CARD 


VARIABLES 


FORMAT 


I 


i 


1 LABEL  8A10 

2 NCHAR,  N,  LAMBDA,  315,  2F10.5 

DU,  L 

3 QMIN,  QMAX,  DQ1,  7F10.2 

ETAMIN , ETAMAX, 

DETA1 , DETA2 

Values  of  user  input  data  used  to  complete  a successful  run  of  the 
BUBPLT  program  on  the  CDC  6000  Series  Computer  are  provided  below: 

LABEL  = FROM  MODELS  1094CV  AND  1094  - B2  AT  V=5.360  FT/SEC  (AUG.  1968) 

NCHAR  =72 

N = 100 

LAMBDA  = 0 

DU  =0.1 

L = 14.01 

QMIN  = - . 08 

QMAX  = . 08 

DQ1  = .01 

ETAMIN  = .60 

I 

ETAMAX  =2.0 
DETA1  = .05 
DETA2  = . 20 

Table  9 provides  a listing  of  this  sample  input  data  in  the  format 
required  by  program  BUBPLT. 

The  input  variables  FH  and  GH  are  read  from  device  number  3 and  should 
be  the  free-wave  spectra  output  of  program  WAVECT  without  a bulb.  Input 
variables  FBO  and  GBO  are  read  from  device  number  4 and  should  be  the  free- 
wave  spectra  output  of  program  WAVECT  with  a bulb. 


Table  9 - Sample  Input  for  Program  BUBPLT 


1 

1 


3 


The  BUBPLT  program  has  printed  output  consisting  of  the  maximum  and 
minimum  of  the  P and  Q values,  as  well  as,  eta  plus  plotter  commands  output 
to  a magnetic  tape  which  will  drive  the  CalComp  plotter.  A sample  of  the 
CalComp  plots  for  program  BUBPLT  is  shown  in  Figure  4. 

The  input  and  output  of  program  BUBPLT  is  performed  via  devices  number 
3,  4,  5,  6,  and  7.  Device  number  5 should  be  assigned  to  the  card  reader 
for  the  user  inputs,  device  6 to  the  line  printer  for  the  printed  output, 
device  7 to  a 7-track  magnetic  tape  for  the  CalComp  plotter  commands,  and 
device  numbers  3 and  4 should  be  assigned,  respectively,  to  the  permanent 
file  containing  the  free-wave  spectra  output  of  program  WAVECT  run  with 
and  without  a bulb. 

Sample  control  cards  to  complete  a successful  run  of  the  BUBPLT 
program  on  the  CDC  6000  Series  Computers  are  listed  below: 

CHREPLOT , CM60000 , T120 ,MT1 , P2 . 

CHARGE , CHRE , XXXXXXXXXX . 

FTN(T) 

ATTACH , TAPE3 , PERMFILE3 , ID=CHRE . 

ATTACH , TAPF.4 , PERMFILE4 , ID=CHRE . 

VSN,TAPE7=SL0T33=RAE03 . 

REQUEST, TAPE7, HI, RING.  (SL0T33.RAE03) 

ATTACH, CALC9 36. 

ATTACH, CALCFN. 

LDSET , LIB=CALC936 . 

LDSET , LIB=CALCFN . 

LGO. 

RETURN, TAPE 7. 

7/8/9  END  OF  RECORD  CARD 

Source  Deck 

7/8/9  END  OF  RECORD  CARD 
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Data  Cards 


6/7/ 8/9  END  OF  FILE  CARD 

Note:  PERMFILE3  AND  PERMFILE4  are,  respectively,  the  names  of  the  permanent 
files  containing  the  free-wave  spectra  output  of  program  WAVECT  with 
and  without  a bulb. 

COMPUTER  PROGRAM  DOCUMENTATION 

This  section  of  the  report  contains  the  calling  sequences  and  external 
references  of  the  five  computer  programs  and  related  subroutines  which 
collectively  will  be  used  for  the  analysis  of  longitudinal  wave  cuts.  The 
documentation  in  this  section  describes  the  programs  as  used  on  the  CDC  6000 
series  computers  at  DTNSRDC.  This  section  is  designed  to  assist  in  the 
program  maintenance  function. 

A summary  of  the  external  references  is  outlined  below.  Following 
the  summary,  documentation  for  programs  and  non-system  functions  and  sub- 
routines is  provided. 


NONDIM  (no  references) 


NAME: 

PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 


SUBROUTINES  CALLED: 
COMMENTS : 


PROGRAM  NONDIM 

To  convert  wave  cut  information  from  analog  units 
to  feet,  and  then  nondimensionalize  the  information 
for  the  wave  analysis  program 

MAIN  PROGRAM 

TAPE5  - Card  Reader 

TAPE6  - Line  Printer 

TAPE4  - Disk  File  (Nondimensional  Data) 

None 

None 


- 


NAME: 

PURPOSE: 


CALLING  SEQUENCE: 
ARGUMENTS : 


I 


( i 


SUBROUTINES  CALLED: 
FUNCTIONS  CALLED: 
COMMENTS : 


r i 


PROGRAM  WAVECT 

To  perform  wave  analysis  using  the  fourier 
transform  method  with  an  optional  truncation 
method  which  can  be  applied  to  height  data 

MAIN  PROGRAM 

TAPES  - Card  Reader 

TAPE 6 - Line  Printer 

TAPE3  - Disk  File  (Free  Wave  Spectra  Without  Bulb) 
TAPE4  - Disk  File  (Nondimensional  Data) 

CONST,  TRANS,  TRUNC,  FRES,  FANDG,  RESIST 

SQRT 

None 


! 
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NAME: 


SUBROUTINE  CONST 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 

j 

I 


I 


■ 

I 

j / 

SUBROUTINES  CALLED: 
FUNCTIONS  CALLED: 
COMMENTS : 

I 


i 


By  using  least  square  analysis,  CONST  determines 

the  asymptotic  behavior  of  the  waves  behind  the 

ship:  ZETA  - (Cl*Cos  (X)  + C2  *Sin  (X)  )/  (C3  - X) , 

as  X approaches  (-)  infinity. 

CALL  CONST  (ZETA,  M,  MP,  XONE,  DX,  Cl,  C2,  C3) 

ZETA  - The  vector  of  length  MAX  containing  the  input 
points  for  the  longitudinal  wave  profile 

M - An  integer  Indicating  the  number  of  points 

which  are  to  be  analyzed  from  the  wave  profile 

MP  - An  integer  which  indicates  at  which  point  on 
the  wave  profile  the  analysis  for  determining 
the  asymptotic  behavior  of  the  wave  profile 
should  begin  - MP  is  not  read  unless  MU  = 0. 

XONE  - The  distance  of  the  first  input  point  forward 
of  the  origin  on  the  model 

DX  - The  distance  between  the  input  in  the  wave 
profile 

Cl  - The  real  part  of  the  complex  amplitude  of  the 
asymptotic  wave  behavior 

C2  - The  imaginary  part  of  the  complex  amplitude 
of  the  asymptotic  wave  behavior 

C3  - A nondimens ional  distance  indicating  where  the 
asymptotic  behavior  of  the  wave  begin  in 
relation  to  the  origin  on  the  model  - C3  is 
read  only  if  MU  = 0. 

None. 

COS,  SIN,  SQRT 


NO"“ 


NAME: 


SUBROUTINE  TRANS 


PURPOSE: 

CALLING  SEQUENCE: 

ARGUMENTS : 


SUBROUTINES  CALLED: 


To  compute  the  sine  and  cosine  components  of  the 

Fourier  Transform  at  transverse  wave-number  SNU 

CALL  TRANS  (ZETA,  MU,  M,  DX,  XM,  SNU,  TNU,  CNUSTR, 

SNUSTR) 

ZETA  - The  vector  of  length  MAX  containing  the 
input  points  for  the  longitudinal  wave 
profile 

MU  - An  Integer  control  character  indicating 
whether  slope  or  height  data  is  to  be 
read  and  whether  or  not  the  truncation 
correction  is  to  be  made  in  the  case  of 
height  data 

MU  < 0,  height  data  without  truncation 
correction 

MU  * 0,  height  data  with  truncation 
correction 

MU  >0,  slope  data 

M - An  integer  indicating  the  number  of  points 

which  are  to  be  analyzed  in  the  wave  profile 

DX  - The  distance  between  the  input  points  in  the 
wave  profile 

XM  - The  point  at  which  the  asymptotic  behavior 
of  the  wave  is  assumed  to  start  - point 
from  which  the  truncation  correction  is 
made 

SNU  - The  longitudinal  wave  number  - Sec  0 

TNU  - Tan  0 

CNUSTR  - The  cosine  component  of  the  Fourier  Transform 
of  the  wave  profile  without  the  truncation 
correction 

SNUSTR  - The  sine  component  of  the  Fourier  Transform 
applied  to  the  wave  profile  without  the 
truncation  correction 

None 
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NAME: 


SUBROUTINE  TRUNC 


I 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 


To  determine  the  truncation  correction  to  the  free- 

wave  spectra  due  to  finite  length  of  record 

CALL  TRUNC  (XM,  SNU,  Cl,  C2,  C3,  DCSTR,  DSSTR) 

XM  - The  point  at  which  the  asymptotic  behavior 
of  the  wave  is  assumed  to  start  - point 
from  which  the  truncation  correction  is  made 

SNU  - The  longitudinal  wave  number  - Sec  0 

Cl  - The  real  part  of  the  complex  amplitude  of  the 

asymptotic  wave  behavior 

C2  - The  imaginary  part  of  the  complex  amplitude 
of  the  asymptotic  wave  behavior 

C3  - A nondimensional  distance  indicating  where 
the  origin  for  use  determining  where  the 
asymptotic  behavior  of  the  wave  begin  in 
relation  to  the  origin  on  the  model 

DCSTR  - The  truncation  correction  for  the  cosine 
component  of  the  Fourier  Transform  of  the 
wave  profile  without  the  truncation  correction 

DSSTR  - The  truncation  correction  of  the  sine 

component  of  the  Fourier  Transform  which 
applies  to  the  wave  profile  without  the 
truncation  correction 


SUBROUTINES  CALLED: 

FRES 

FUNCTIONS  CALLED: 

COS, 

COMMENTS : 

None 

SIN,  SQRT 


i 


i 


r * 
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j 


i 


NAME: 

SUBROUTINES  FRES 

PURPOSE: 

To  calculate  the  value  of  the  Fresne  integral  from 
zero  to  Z and  return  the  real  (CT)  and  imaginary  (£>T) 
parts  using  the  method  developed  by  Boersma  (1960)^ 

CALLING  SEQUENCE: 

CALL  FRES  (Z,  CT,  ST) 

ARGUMENTS : 

Z - The  upper  limit  of  the  Fresnel  integral 

CT  - The  real  part  of  the  Fresnel  integral  from 
zero  to  Z 

ST  - The  imaginary  part  of  the  Fresnel  integral 
from  zero  to  Z 

SUBROUTINES  CALLED: 

None 

FUNCTIONS  CALLED: 

CEXP , CSQRT 

COMMENTS : 

None 

/ 


NAME: 


SUBROUTINE  FANDG 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 


SUBROUTINES  CALLED: 
FUNCTIONS  CALLED: 


To  compute  the  "true"  free-wave  spectra,  F and  G. 

CALL  FANDG  (MU,  Y,  UNU,  VNU,  CNUSTR,  SNUSTR,  F,  G) 

MU  is  an  integer  control  character 

MU  < 0;  Height  data  without  truncation  correction 
MU  = 0;  Height  data  with  truncation  correction 
MU  > 0;  Slope  data 

Y - The  distance  from  the  center  line  of  the 

model  to  the  center  of  the  wave  wire 

UNU  - The  transverse  wave  number  - Sec  8 • Tan  0 

VNU  - The  function  of  wave  number  used  in 

evaluating  various  integrands  ~ VNU  = 

^(1  + 4 Sec  0 Tan  6) 

CNUSTR  - The  cosine  component  of  the  Fourier  transform 
of  the  wave  profile  without  the  truncation 
correction 

SNUSTR  - The  sine  component  of  the  Fourier  transform 
applied  to  the  wave  profile  without  the 
truncation  correction 

F - The  sine  component  of  the  free-wave  spectra 

G - The  cosine  component  of  the  free-wave  spectra 

None 

COS,  SIN 


COMMENTS : 


None 


NAME: 


SUBROUTINE  RESIST 


m 


! 


>>  - 


PURPOSE: 

CALLING  SEQUENCE: 

ARGUMENTS : 


To  compute  the  incremental  increase  in  wave  resistance 
and  side  force  due  to  a single  step  in  wave  number 

CALL  RESIST  (DU,  VNU,  SNU,  UNU,  CNUSTR,  SNUSTR,  DR, 

DT) 

DU  - The  size  of  the  increment  for  the  transverse 
wave  number  which  is  used  to  determine  the 
N elementary  waves  which  are  to  be  analyzed 

VNU  - The  function  of  wave  number  used  in  evaluating 
various  integrands.  VNU  - ^(1  + 4 Sec0  Tan0) 

SNU  - The  longitudinal  wave  number  - Sec0 

UNU  - The  transverse  wave  number  - Sec9  • Tan0 


CNUSTR  - The  cosine  component  of  the  Fourier  transform 
of  the  wave  profile  without  the  truncation 
correction 

SNUSTR  - The  sine  component  of  the  Fourier  transform 
applied  to  the  wave  profile  without  the 
truncation  correction 


DR 

- The  incremental 
due  to  a single 

increase  in  wave  resistance 
step  in  wave  number 

DT 

- The  incremental 
to  a single  step 

increase 
in  wave 

in  side  force  due 
number 

SUBROUTINES  CALLED: 

None 

COMMENTS : 

None 

49 


1.  ' 


m 


NAME: 

PURPOSE: 


CALLING  SEQUENCE: 
ARGUMENTS : 


SUBROUTINES  CALLED: 
FUNCTIONS  CALLED: 
COMMENTS : 


PROGRAM  SPCTRA 


To  plot  the  sine  and  cosine  component  of  the  free- 
wave  spectra 


MAIN  PROGRAM 

TAPE4  - Disk  File  (Free-Wave  Spectra) 

TAPE5  - Card  Reader 

TAPE6  - Line  Printer 

TAPE7  - 7-Track  Magnetic  Tape 

NEWMIN,  AXIS,  PLOTS,  LINE,  PLOT,  SCALE,  SYMBOL 

SQRT 

None 


NAME: 

PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 


SUBROUTINES  CALLED: 
COMMENTS : 


SUBROUTINE  NEWMIN 

To  determine  the  scaling  to  be  used  when  plotting 
the  free-wave  spectra 

CALL  NEWMT“  (AMIN,  AX,  DV,  AORG) 

AMIN  - The  vertical  range  in  which  the  spectra  is 
to  be  plotted 

AX  - The  maximum  value  of  the  spectra  arrays  F, 
G,  and  E 

DV  - The  number  of  data  units  per  inch  of  axis 
AORG  - A shifting  factor  for  the  spectra  arrays 
None 
None 
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A 


NAME: 


PROGRAM  BUBOPT  (PRINTER) 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 


SUBROUTINES  CALLED: 
FUNCTIONS  CALLED: 
COMMENTS : 


To  determine  coarse  contours  of  the  bulb  influence 
factor 

MAIN  PROGRAM 

TAPE3  - Disk  File  (Free-Wave  Spectra  Without  Bulb) 

TAPE4  - Disk  File  (Free-Wave  Spectra  With  Bulb) 

TAPE5  - Card  Reader 

TAPE6  - Line  Printer 

CONTOR 

SUM,  SIN,  COS,  SQRT 
None 
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NAME: 


SUBROUTINE  CONTOR 


PURPOSE: 

CALLING  SEQUENCE: 

ARGUMENTS : 


SUBROUTINES  CALLED: 
FUNCTIONS  CALLED: 


To  determine  the  contours  of  constant  ETA  in  the 
P-Q  space  and  print  them 

CALL  CONTOR  (QMIN,  QMAX,  DQ,  ETAMIN , ETAMAX,  DETA, 

A,  B,  C,  I) 

QMIN  - The  starting  value  of  Q 
QMAX  - The  ending  value  of  Q 

DQ  - The  transverse  size  for  Q 

ETAMIN  - The  minimum  ETA  for  which  contours  are  wanted 

ETAMAX  - The  maximum  ETA  for  which  contours  are  wanted 

DETA  - The  ETA  distance  between  contours 

A,B,C  - For  a given  Q,  ETA  is  quadratic  in  P 

I - A counter  used  to  dimension  the  argument 

B in  the  subroutine  CONTOR 

None 

SQRT 


COMMENTS : 


None 


If 


r 


NAME: 

PURPOSE : 

CALLING  SEQUENCE: 
ARGUMENTS : 


SUBROUTINES  CALLED: 
COMMENTS : 


FUNCTION  SUM 

To  determine  the  interference  wave  resistance  of  a 
hull  bulb  combination  - This  is  used  to  evaluate  the 
B term  in  the  quadratic  expression  (in  bulb  volume) 
for  wave  resistance  at  a given  longitudinal  position. 

S - SUM  (NN,  DU,  FA,  FB,  GA,  GB,  W) 

NN  - The  number  of  free-wave  spectra 

DU  - The  transverse  wave  number  step  size 

FA  - The  sine  component  of  the  free-wave  spectra 
for  the  hull 

FB  - The  sine  component  of  the  free-wave  spectra 
for  the  bulb 

GA  - The  cosine  component  of  the  free-wave  spectra 
for  the  hull 

GB  - The  cosine  component  of  the  free-wave  spectra 
for  the  bulb 

W - Denominator  of  the  Integrand 

W - + 4U2)  / (1  + V 1 + 4U2) 

None 

None 


NAME: 

PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 


SUBROUTINES  CALLED: 
FUNCTIONS  CALLED: 
COMMENTS : 


PROGRAM  BUBPLT 

To  determine  and  plot  the  contours  of  the  constant 
bulb  influence  factor  (ETA) 

MAIN  PROGRAM 

TAPE3  - Disk  File  (Free-Wave  Spectra  Without  Bulb) 

TAPE4  - Disk  File  (Free-Wave  Spectra  With  Bulb) 

TAPE5  - Card  Reader 

TAPE6  - Line  Printer 

TAPE7  - 7-Track  Magnetic  Tape 

VOL,  PLTST , PLOT 

SUM,  INFLEX,  XINTRP , COS,  SIN,  SQRT 
None 
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NAME: 

PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 


SUBROUTINES  CALLED: 
FUNCTIONS  CALLED: 
COMMENTS: 


SUBROUTINE  VOL 

2 

To  solve  the  quadratic  equation,  ETA  = AP  + BP  + C, 
for  the  bulb  volumes  which  yield  a given  value  of  ETA 

CALL  VOL  (ETA,  A,  B,  C,  PA,  PB,  T) 

ETA  - The  given  value  of  bulb  influence  factor  for 
which  contours  are  designed 

A - Coefficient  on  quadratic  term  in  bulb  volume 

B - Coefficient  on  linear  term  in  bulb  volume 

C - Coefficient  on  constant  term 

PA  - The  lower  bulb  volume  for  the  contour 

PB  - The  upper  bulb  volume  for  the  contour 

T - Used  as  a switch  on  the  number  of  values  of  P 

None 

SQRT 

None 
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NAME: 


SUBROUTINE  PLTST 


I 


I 


I 


i 


PURPOSE:  To  perform  the  actual  calcomp  plotting  of  the  contours 

of  the  constant  bulb  influence  factor 

CALLING  SEQUENCE:  CALL  PLTST  (N,  A,  AX,  B,  BX,  BMIN,  K,  L,  ETA,  LABEL, 

NCHAR,  IENTRY) 

N - The  number  of  points  on  a given  contour  to 

be  plotted 

A - The  Q values  to  plot  on  a contour 

AX  - The  scaling  constant  for  the  Q plot 

B - The  P values  corresponding  to  Q 

BX  - The  scaling  constant  for  the  P plot 

BMIN  - The  minimum  value  of  P 

K - A switch  used  to  denote  open  or  closed  contours 

L - A switch  used  to  denote  solid  or  dashed  lines 

ETA  - The  value  of  the  bulb  influence  factor  being 
plotted 

LABEL  - A title  for  the  output 

NCHAR  - The  number  of  characters  in  the  title 

IENTRY  - A switch  used  to  denote  scale  set-up,  plot 
contours  or  dashed  lines 

SUBROUTINES  CALLED:  PLOTS,  PLOT,  SYMBOL,  AXIS,  LINE,  NUMBER,  DASHL,  LINE1 

FUNCTIONS  CALLED:  INFLEX 

COMMENTS : None 


NAME: 


FUNCTION  SUM 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS: 


SUBROUTINES  CALLED: 
COMMENTS : 


To  determine  the  interference  wave  resistance 
of  a hull  bulb  combination  - This  is  used  to  evaluate 
the  B term  in  the  quadratic  expression  (in  bulb 
volume)  for  wave  resistance  at  a given  longitudinal 
position. 

S - SUM  (NN,  DU,  FA,  FB,  GA,  GB,  W) 

NN  - The  number  of  free-wave  spectra 

DU  - The  transverse  wave  number  step  size 

FA  - The  sine  component  of  the  free-wave  spectra 
for  the  hull 

FB  - The  sine  component  of  the  free-wave  spectra 
for  the  bulb 

GA  - The  cosine  component  of  the  free-wave  spectra 
for  the  hull 

GB  - The  cosine  component  of  the  free-wave  spectra 
for  the  bulb 

W - Denominator  of  the  Integrand: 

W « V 1 + 4U2  / 1 W 1 + 4U2 


None 

None 
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NAME 


FUNCTION  INFLEX 


n 


' 


! 

I 


PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 

SUBROUTINES  CALLED: 
COMMENTS: 


To  determine  the  inflection  point  in  a given  set  of 
data  - the  vector  A 

I - INFLEX  (A,  N) 

A - A vector  of  values  which  has  a maximum  or  minimum 

N - The  number  of  elements  in  the  vector  A 

None 

None 
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.JPC 


NAME: 

PURPOSE: 

CALLING  SEQUENCE: 
ARGUMENTS : 


SUBROUTINES  CALLED: 
COMMENTS : 


FUNCTION  XINTRP 

To  interpolate  a set  of  data  Y for  a value  of  Z = X 
using  quadratic  interpolation 

XX  - XINTRP  (X,  Y,  N,  Z) 

X - The  vector  of  X values  in  which  to  interpolate 

Y - The  vector  of  Y values  in  which  to  interpolate 

N - The  number  of  elements  in  vectors  X and  Y 

Z - The  value  of  X at  which  a value  of  Y is  desired 

None 

A value  of  100.0  is  returned  to  denote  interpolation 
fails 


r 


APPENDIX  A - NONDIMENSIONAL  SCHEME 


In  order  that  the  results  of  the  Wave  Analysis  would  be  independent 
of  the  size  of  the  model  or  ship  from  which  the  data  was  obtained,  it  was 
decided  that  nondimensional  quantities  would  be  used  in  the  actual  wave 
analysis.  For  use  in  the  dimensional  analysis,  a set  of  basic  units  of 
density  (p) , gravity  (4),  and  velocity  (V)  were  chosen,  for  use  in  repre- 
senting all  data  in  a nondimensional  form. 

We  can  now  form  nondimensional  quantities  as  follows;  let  Q be  any 
dimensional  quantity  such  that: 

[fi]  - [MaLbTC] 

Then,  letting  Q be  the  corresponding  nondimensional  quantity,  we  have: 


a 


With  tie  following  relationships  between  V,  and  M,  L,  T: 


General 

(English) 

(Metric) 

[V] 

L/T 

(ft/sec) 

(m/s) 

[£.] 

L/T2 

(ft/sec2) 

(m/ s2 ) 

t£l 

M/L3 

(slug/ft3  = Lb  sec2/ft4) 

(Kg/m5  ) 

we  can  obtain  the  following  equation  relating  a,  b,  c and  a,  8,  Y 
MaLbTC  = (ML~3)a(LT"2)e(LT_1)Y 
which  yields  the  following  relationships: 


b “ -3a  + 8 + Y 
r « -2ft  - Y 


APPENDIX  B - MATHEMATICAL  DEVELOPMENT* 


According  to  theoretical  hydrodynamics  the  potential  for  a source 
distribution  over  some  region  (D)  in  the  lower  half  plane  (z  < 0)  can  be 
described  by  the  following  integral 


<P(x,y,z) 


= j o(x' ,y' ,z')G(x,y,z,x' ,y' ,z')dD  . 
D 


(1) 


where  a(x',y',z')  is  the  local  source  strength  for  a point  (x',y',z'); 
and  the  Green's  function  G(x,y,z,x' ,y' ,z')  is  the  potential  of  a source 
of  unit  strength  located  at  a point  with  coordinates  (x1,  y',  z'). 
Havelock  (1932)6  gives  the  following  expression  for  the  Green's  function: 


G(x,y,z,x' ,y' ,z') 


r - r 
1 2 

r r 
1 2 


+ lim 
JZ-’O 


00 

i/-/ 

-ir  n 


dk 


sec2  9 exp{k(z  + z'  + iui)} 
k - sec2  0 + iy  sec  0 


(2) 


where  r^  = [(x  - x')2  + (y  - y')2  + (z  - z')2]  , 

r = [(x  - x')2  + (y  - y')2  + (z  + z')2]  , and  (3) 

2 

co  = (x  - x')  cos  0 + (y  - y')  sin  0 . 

Defining  the  steady  state  free  surface  deformation  as; 

z = £(x,y)  , 

we  can,  by  the  assumption  of  small  amplitude,  define  this  free  surface 
deformation  in  terms  of  the  potential  function  as: 

z = £(x,y)  = <(>x(x,y,0)  (A) 

If  we  define  the  Kochin  functions  H(k,0)  and  J(u,s)  as  follows: 

H(k,0)  = 4tt J"  a(x'  ,y'  ,z')exp[kz'  + ik(x  cos  0 + y'  sin  9)]dD  and  (5) 
D 

*For  a more  detailed  analysis  see  Sharma  (1966) 
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(u,s)  - 4ir  a(x' ,y' .z'Jexpfz'  ^ uz  + sz  + i(sx'  + uy')]dD  . (6) 


Then  we  can  obtain  the  following  relationships  between  H and  J: 


H(k,0)  = J(k  sin  8,  k cos  8)  and 


J(u,s) 


« (V”u 


2 + s2,  tan 


Using  the  Kochin  functions,  we  can  now  obtain  the  free  surface 


deformation  as: 


C(x,y)  = Re 


} f 

lim  — / d8  / c 

H-0  4tt2  J J 


H(k,8)k  sec  8 


k - sec2  0 + ip  sec  0 


exp[ik(x  cos  8 + y sin  9)]  , 


where  H(k,0)  is  the  complex  conjugate  of  H(k,8) 


Defining  the  wave  resistance  Rw  by  the  following  equation. 


Rw  * +4',r  / a(x,y,z)4>x(x,y,z)dD  , 
D 


and  substituting  for  4>x(x»y»z)»  we  obtain 


w “ 2tt  / 


H(sec2  0,0)  ‘H(sec2  0,0)sec3  0 d 0 . 


From  the  above  equation  for  the  wave  resistance  it  can  be  observed 
that  the  Kochin  function  has  degenerated  from  a function  of  two  variables 
(k,0)  to  what  is  essentially  a function  of  one  variable  (0).  From  this 


it  can  be  seen  that  k must  be  defined  as  follows: 


k = 8ecz  0 . 


Outside  the  region  near  the  source  distribution,  the  surface 
deformation  Is  caused  predominantly  by  those  waves  for  which  k = secs  8. 
These  are  the  so-c#lled  free-waves. 

It  can  be  shown  that,  for  any  y,  as  x -*■  00  the  surface  deformation 
becomes : 


«*,y>  f - 

2 i VI  + 4u2 


(1  + \l  + 4us ) du  , 


and  that,  for  any  x,  as  Y -*•  00 


C(x,y)  = I 


1 f J (—is  VI 

* J ~rr 


- s ,s) 


1 - s' 


exp[lsx  ■ s Vl  • ssy]3s  ds 


if 


J(s\s2  - l,s) 


's2  - 1 


exp[i(sx  + sys2  - l,y)]s2ds  , 


where  s is  related  to  u by  the  following  relationship: 


I 1/2 


From  this  it  can  be  seen  that  the  various  forms  of  the  Kochin  function 
are  related  as  follows: 


H(sec2  9,0)  = J(u,s)  = J(sys 2 - 1,  |s|)  . I 

From  this,  it  can  be  seen  that  the  variables  u,  s,  and  0 are  related  as 
follows : 

u = Sec  9 Tan  8 , 
s = Sec  9 , and 

u = s Vs3  - l . ( 
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Using  these  variables  we  can  define  a "free-wave  spectra"  in  terms  of 
the  two  real  functions,  F and  G,  of  the  variable  u, 


G(u)  + iF(u)  « 2(1  1=^L-4u— ) J(u,s)  , 

VI  + 4u2 


and  the  wave  elevation  as  x -*■  00  can  be  found  by  the  Fourier  integral. 


C(x,y) 


l.  r 

4tt  J 


[F(u)  sin  (sx  + uy)  + G(u)  cos  (sx  + uy)]du 


Taking  a longitudinal  wave  cut  (y  held  constant),  the  following 
Fourier  transform  can  be  applied  to  the  wave  profile: 


C(s,y)  + iS(s,y)  = J £(x,y)eisx  dx  . 


If  the  value  of  y has  been  chosen  such  that  it  is  great  enough  that 
the  waves  are  free  waves,  we  can  apply  the  Fourier  inversion  theorem  to 
equation  (19),  and  for  1 < s < 00  obtain  the  following  relationships: 

Re[J(sVs2  - l,s)exp{is>/ s2  - 1 y}]  * — i C(s,y)  and 

s2 

I [J(s Vs2  - l,s)exp{-isVs2  - 1 y}]  ^ S(s,y)  . 

m s2 

From  these  equations  we  can  easily  solve  for  the  free-wave  spectra, 
in  terms  of  the  Fourier  transform  of  equation  (19): 

F(u)  ■ [C(s,y)  sin  (uy)  + S(s,y)  cos  (uy)  ] and 

(2s2  - 1 

G(u)  - ^ 8 ~ — [C(s,y)  cos  (uy)  + S(s,y)  sin  (uy) ] • 


These  equations  (21)  are  true  for  height  data;  however,  there  can 
be  significant  errors  from  using  these  equations  with  only  a finite 
length  of  wave  cut.  To  avoid  these  errors,  it  is  possible  to  instead 
measure  the  transverse  wave  slope  h^  which  would  yield  the  following 
equations  for  the  free-wave  spectra. 


4 

F(u)  = [C  (s,y)  cos  (uy)  - S (s,y)  sin  (uy)  ] 

s (2s2  - 1)  y y 

4 

G(u)  = [C  (s,y)  sin  (uy)  + S (s,y)  cos  (uy)] 

s(2s2  - 1)  y y 


(22) 


A method  which  might  possibly  be  used  to  correct  for  the  use  of  a 
finite  length  of  longitudinal  cut,  with  height  data,  would  be  to  assume 
that  at  a large  distance  from  the  origin,  the  waves  decay  asymptotically, 
and  that  they  may  be  fit  by  an  equation  such  as  follows,  for  x -*•  -®: 


C(x,y) 


(C  cos  x + C sin  x) 

l 2 

■J  C - x 

’ 3 


(23) 


where  the  unknown  constants  C , C , and  C might  be  determined  by  a 

x 2 3 

least-squares  fit. 


In  order  to  determine  the  error  due  to  truncation  let's  define  a 
modified  Fourier  transform, 

00 

C*(s,y)  + iS*(s,y)  - J * V s2  - 1 C(s,y)elsx  dx  . (24) 

—00 


Then,  if  we  define  the  point  x^  as  the  point  at  which  the  wave  profile 
has  been  truncated,  the  error  becomes: 


AC*(x,y)  + iAS*(s,y) 


2 M i 

= J V s2  - 1 C(x,y)e  sx  dx  , 


(25) 


where  the  values  of  C(x,y)  are  defined  by  using  equation  (23).  Using 
equations  (25)  and  (23)  we  can  obtain  the  values  of  AC*  and  AS*  by 
using  the  following  equations: 
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where 


AC*  * [d^Cp(z  } + d^Sp(z  ; + d^Cp(z”)  + d^z")]  and 
AS*  " VF td2CF(z+)  " d1SF(z+)  + d„CF(z)  " d3SF(z_)]  ’ 


d 

1 

d 

2 
d 

3 
d 

4 


^ ir-  1[C-  cos  (s  + 1)  + C 

Vtt-  1[C  sin  C (s  + 1)  - C 

1 3 

•JirT"T[C  cos  C (s  - 1)  - C 


1 3 

TT+  1[C  sin  C (s  - 1)  + C 


sin  C (s  + 1) ] > 

3 

cos  C (s  + 1) ] , 

3 

sin  C (s  - 1)]  , and 
3 

cos  C (s  - 1)]  . 


C and  S_  are  the  real  and  imaginary  parts  of  the  Fresnel  integral, 
t F 


00  i 1 t2 

Cp(z)  + i Sp(z)  «/  e 2 dt  , 


and 


| 2(C  - Xj^)  (s  + 1)  j»/2 


Using  the  real  and  imaginary  components  of  the  modified  Fourier 
transform  (24)  we  can  now  define  the  wave  resistance  and  transverse 
forces  as 

R - i f <C*>  4 d»  .„d 


APPENDIX  C - NUMERICAL  ANALYSIS  AND  THE 

DEVELOPMENT  OF  THE  ALGORITHM 

We  will  use  the  following  definitions  in  our  evaluation  of  the 
equations  in  Appendix  B: 

V * 0,1,2,...N  the  number  of  an  elementary  wave  with  angle  0 , 


Uy  - V ' Au 

Vv  - (1  + Au^1/2 

Sv  " {(1  + V/2}1/2 


Tv  " uv/sv 
Ej  " 1/2 


( sec  9V  tan  0y  ) , 
(1  + 2 tan2  0y  ) , 
(sec  0V)  , 

(tan  0V)  , 


j = j . and  j 
J Jmin  Jmax 

1 - 1 . + l,...,j  - 1 

J Jmin  ,Jmax 


Ax  the  spacing  between  points  in  the  wave  profile,  and 
M = the  number  of  the  most  negative  x-location  in  the  wave  profile. 
Taking  the  truncated  Fourier  transform  integral: 


C*(s)  + iS*(s)  = VSv  « , 


we  can  separate  the  real  and  imaginary  parts  and  integrate  by  the  trapi- 
zoidal  rule  obtaining  the  equations: 


SV*  " TvlAxl^  sin  (XjSv)ej  1 

cv* = TviAxi£  s c°s  (xjVej  • 


Using  these  results  with  equation  (21)  we  can  obtain  the  free-wave 
spectra: 
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(35a) 


I 


Fv  “ 7~{V  8ln  (UvY)  + Sv*  cos  (UvY)}  ’ 


Gv  ■ v“^cv*  cos  <UVY)  ■ Sv*  Sin  <V>  > » and 


(35b) 


Ev  " {<V2  + <Gv)2}1/2 


(35c) 


Then  by  integrating  equations  (30)  and  (31)  by  the  trapizoidal  rule  we 
can  obtain  the  wave  resistance  and  transverse  force: 


3s  L + 

1=0  J 


■ v 


V 


and 


(36) 


J-o 


*)!  + <V,M 


(37) 


When  being  used  with  a slope  cut  the  equations  develop  as  follows: 


C*(s)  + i S*(s) 


oo 

kf 


Cy(x»y)e1SXdx  . 


(38) 


m 


This  equation  then  separates  into  the  sums; 


V 


|Ax| 
s. 


m 


FE  Cy  sin  (XjSv)£j  and 

V J-l  J 


(39) 


JM£ 


V 


2^  S,  COS  (xjsv)e3  • 

v j«l  j 


The  free  wave  spectra  are  given  by  the  following  equations: 
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A 


Fv  “ v~  {Cv*  cos  (UvY)  - sv*  8in  (UVY)}  and 


Gv  = T-  {cv*  sin  (UvY)  + Sv*  cos  (UvY)}  • 

V 


The  values  of  Ey,  R^,  and  are  then  calculated  as  in  equations 
(35c),  (36),  and  (37)  respectively. 

The  evaluation  of  Cv*  and  Sy*  can  be  much  simplified  if  it  is  noted 
that  the  summations  of  equations  (34)  and  (39)  are  simply  Chebyshev 
polynomials  to  which  the  following  recursive  formulas  apply: 


UM+1  UM+2  = 0 


where 


Un  = ?n  + 2 COS  (SvAx)Un+l  " Un+2  * 


C = C /2  for  n = 0,M 
n n 


This  equation  is  iterated  from  M to  zero,  and  S^*  and  C^*  are  calculated 
from  the  following  equations: 


Sv*  = TjAX|U  sin  (SyAX)  and 
V = TV|AX1^U°— -U-2  ) . 


The  truncation  corrections  ACy*  and  AS^*  are  calculated  and  added 

to  Cy*  and  S^*  from  equation  (34)  to  obtain  C^*  and  S^*,  and  then  the 

free-wave  spectra,  transverse  force  and  wave  resistance  are  calculated 

using  equations  (35),  (36),  and  (37).  The  truncation  corrections  ACy* 

ASy*  are  calculated  using  equations  (26)  and  (27),  where  the  values  of 

the  constants  C , C , and  C are  determined  as  follows:  C is  chosen 
12  3 3 

as  the  nondimensional  distance  from  the  origin  on  the  model  to  the  origin 
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of  the  most  significant  waves  on  the  hull.  The  constants  Cj  and  C2 
are  calculated  by  least  squares. 


The  error  6j  is  calculated  such  that 


C cos  x.  + C sin  x. 
Kmc-  -1 j 2 1 

3 3 >/vrV 


Then  the  sum  of  the  squared  errors,  Q,  is  calculated. 


£ <va  • 


j-M 


It  is  desired  to  find  C and  C such  that  Q is  a minimum;  and 

1 2 

Q is  a minimum  when 


3C 

1 


M / C cos  X.  + C sin  X.  \ cos  X, 

_V  Ir  - -1 j 2 J-J  — 4-  -0 


and  when 


ac 

2 


M / C cos  X.  + C sin  X.  \ sin  X. 

-v  I c - -1 — . J -J- H . = 4?  - 0 

is  V3  vc  - Xj  - xj 


From  this  we  obtain : 


j-M 


M cos  X. 


C - X, 
f 3 J 


+ c 


£ 


sin  X^  cos 


j-M' 


C - X, 
3 j 


cos  X 


1 


(a) 


(b) 


£■  3 

(c) 


and 
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(43) 


(44) 


(45) 


1 


I 


In  equations  (44),  (45),  and  (46);  M'  corresponds  to  the  point  at  which 
the  wave  profile  is  assumed  to  begin  its  asymptotic  behavior.  Both  M 
and  M*  should  be  chosen  so  that  (Xj^  - X^,)  is  either  one  or  else  one- 
half  of  one  fundamental  wave  length. 

The  Fresnel  integrals  of  equation  (28)  are  evaluated  by  using  the 

4 

method  developed  by  J.  Boersma  (1960). 
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APPENDIX  D - BULB  OPTIMIZATION  SCHEME 


If  we  are  given  a main  hull  spectrum  Fy^  and  and  a hull-bulb 

spectrum  an<j  G^^^poR  v = 0,1*  . . N for  some  definite  bulb 

size  and  location,  along  with  the  transverse  wave  number  step  size  AU,  and 

2 2 

the  nondimensional  hull  "length:  L • / V_  ■ 1/F^  * we  can  Pre<*ict  t^ie 

new  hull-bulb  spectra  (F  ^^AND  q for  various  bulb  sizes  (P)  and 


Longitudinal  positions  (Q) ; and  calculate  a predicted 
influence  factor  (u)  for  these  spectra. 


and  bulb 


If  we  generate  these  values  of  ri  for  a number  of  positions  and  sizes, 
one  can  even  plot  contours  of  efficiency.  These  contours  might  look  as 
shown  in  Figure  D-l: 


T| =1 . 4 


Present  Bulb 


Optimum  Bulb 


Figure  D-l  - Bulb  Contours  of  Efficiency  for  Various 
Bulb  sizes  versus  Longitudinal  Position 
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It  can  also  be  shown  in  a linear  sense  that  these  efficiencies  are 
a quadratic  function  of  P for  a given  Q.  That  is  n = aPa  + bP  + C,  a(Q), 
b - b(Q) , c - c (Q) . You  will  see  later  that  "a"  is  actually  a constant 
depending  on  the  given  hull  and  bulb,  and  that  "c"  is  always  equal  to  one. 

To  develop  the  necessary  equations  one  assumes  that: 


and 


(bo)  _ „ (hbo)  _ F (h) 
V V V 


(bo)  _ r (hbo)  _ _ (h) 


V = 0,  1 , . • . , N 


If  we  shift  the  bulb  spectra  by  a distance  Xg  we  get: 
(bi)  „ _ (bo) 


F *>  F. .'“-ycos(SvXs)  + Gv(bo)sin(SvXs) 


and 


Gv<bl)  ■ Gv<b°>,:os(svx9)  - Fv<b°>si”<svxs)  • 


where  S = [(1  + \l+4lT*)/2] 1/2  and  U = V*AU.  The  effect  of  changing 
the  bulb  size  by  a given  factor  P yields: 


Fv(»  . PFv(M, 


G <b>  - PG  (bl> 
V V 


and 


The  wave  resistance  can  be  calculated  using  the  following  trapezoidal 
integratic  formula: 


W 8lT  Z'n  C1  + Vl  + AU,,2)  ' V V ’ ’ 


with 


v»0  _ ’v 

1/2  for  v-o  and  N 
1 otherwise 


-We  would  use  the  principle  of  linear  superposition  to  obtain  the  and  Gy 
and  follows: 


F - F ^ + PF 

v v v 

G - G ^ + PG 
V V V 


and 
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By  the  analysis  shown,  one  can  find  the  constants  of  the  quadratic 
equation.  The  important  thing  to  note  is  that  in  the  equation  for  "a": 
N 


AU  V 
8ir 


Ev  V1  + 41,1 


R 87r  v=o  (1  + yjl  + 4U2 

w ’ 


we  are  calculating  a sum  which  is  simply  the  wave  resistance  of  the  bulb 
with  a phase  shift  and  because  a simple  phase  shift  will  not  change  the 
wave  resistance,  the  coefficient  "a"  is  a constant. 

Once  the  coefficients  of  the  quadratic  equation  are  known,  one  can 
find  the  value  of  P for  an  optimum  bulb  at  any  Q and  the  line  h=  1. 


P is  the  relative  bulb  size.  That  is,  the  original  bulb  has  size  P = 1 
and  no  bulb  has  size  P = 0.  The  longitudinal  position  Q is  the  decimal 
shift  in  length  of  the  bulb  relative  to  the  length  of  the  hull. 


/ 
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